Driven Transport on Parallel Lanes with Particle Exclusion and Obstruction 
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We investigate a driven two-channel system where particles on different lanes mutually obstruct 
each others motion extending an earlier model by Popkov and Peschel [JJ . This obstruction may oc- 
cur in biological contexts due to steric hinderance where motor proteins carry cargos by 'walking' on 
microtubules. Similarly, the model serves as a description for classical spin transport where charged 
particles with internal states move unidirectionally on a lattice. Three regimes of qualitatively dif- 
ferent behavior are identified depending on the strength of coupling between the lanes. For small 
and large coupling strengths the model can be mapped to a one-channel problem, whereas a new and 
rich phase behavior emerges for intermediate ones. We derive an approximate but quantitatively 
accurate theoretical description in terms of a one-site cluster approximation, and obtain insight 
into the phase behavior through the current-density relations combined with an extremal-current 
principle. Our results are confirmed by stochastic simulations. 



PACS numbers: 05.40.-a,05.70.Ln,87.10.Mn 
I. INTRODUCTION 

Driven diffusive systems are of importance in various 
fields of physics and biology [5J [3J, since they serve as 
simplistic models for biological transport phenomena j3J- 
0, traffic flow [TJ I5HTU]. fast ionic conductors [TT], as 
well as quasi-classical spin-transport [12] . Furthermore 
they provide valuable insights into non-equilibrium sta- 
tistical mechanics. As an example, and in contrast to 
equilibrium systems, their bulk behavior is sensitive to 
the boundaries [13] . Boundary induced phase transitions 
in one dimension may emerge resulting in complex phase 
behavior. 

The most prominent example of driven lattice gases, 
the Totally Asymmetric Exclusion Process (TASEP), was 
originally proposed as a simple model for the motion of 
multiple ribosomes along a mRNA strand during pro- 
tein translation |14] . In this model particles move unidi- 
rectionally along a one-dimensional lattice, provided the 
next site is empty. Exact solutions, e.g. by employing 
the Bethe ansatz [T5] or a matrix product ansatz [TB] , are 
feasible, yet much insight can readily be obtained from 
simple mean-field considerations [3J. 

Intracellular transport constitutes another fascinating 
biological application [17j of driven systems. Here, molec- 
ular motors such as kinesin or dynein, driven by the hy- 
drolysis of adenosine triphosphate (ATP), move unidirec- 
tionally along microtubules [3] . Macromolecules or other 
cellular constituents, which often are too large to diffuse 
fast enough through the crowded cytosol, are carried by 
motor proteins, and are then actively transported to the 
location where they are needed. Recent theoretical stud- 
ies motivated by these processes have investigated the 
influence of attachment and detachment of the motors 



to the microtubules [TBrl2"2"] . extended particles [S3J, the 
influence of defects on the track [24- 26 , and the compe- 
tition between different motor species [27] [28] . Further 
attention has been paid to transport along several cou- 
pled channels where particles move in parallel. This cou- 
pling can be either achieved by allowing lane-switching 
events [12] I29H37] or by a possible influence of a par- 
ticle in one channel on the motion in the other chan- 
nel [T] [33J ■ Here we consider the latter case and investi- 
gate how mutual obstruction of motor proteins on neigh- 
boring lanes, for example stemming from large cargos 
attached to them, affects the transport properties of the 
system. 

Driven diffusive systems may also serve as a description 
for spin transport with possible applications in the field 
of spintronics. For instance, such spin currents flow in a 
chain of quantum dots where electrons are driven by an 
external voltage in a way that only the lowest energy lev- 
els can be occupied [35]. Hence, not more than one elec- 
tron of each internal state is permitted per site and elec- 
trons located in the immediate vicinity repel each other 
due to Coulomb interaction. A model taking Pauli's ex- 
clusion principle into account while ignoring phase co- 
herence has been investigated recently [TH [30] [3T] , yet 
Coulomb blockade has been neglected. In the present ar- 
ticle we focus on the influence of a mutual obstruction 
mimicking for example Coulomb interaction. To iden- 
tify its effect on the collective transport properties in 
the clearest way, we disregard spin-flip events which can 
cause intriguing behavior on their own |12[ I30|, 131] . 

A simple lattice model that incorporates mutual ob- 
struction on two lanes has been investigated by Popkov 
and Peschel pQ. The steric hindrance there is manifested 
in the hopping rates that explicitly depend on the con- 
figuration on the opposing lane. As a consequence of this 
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coupling of the lanes a variety of peculiar phases arises 
which has been neatly rationalized in terms of a cluster 
approximation. Particularly, symmetry breaking, which 
arises even though the boundary conditions are symmet- 
ric, is observed and analyzed. 

In this article we extend the model of Ref. [1] by con- 
sidering asymmetric boundary conditions and rates in- 
stead of reservoirs at the right boundary. We introduce 
the model in Section [TTJ both in the two-lane and in the 
spin-transport picture. In Section |III| we describe the 
stochastic simulations and provide first insights in how 
particle obstruction affects the behavior. Namly, we iden- 
tify three regimes of qualitatively different behavior. In 
Section [TV] we analytically compute the current-density 
relations within a one-site cluster approximation. Sec- 
tion [V] presents a discussion on how the current-density 
relation, obtained via an extremal-current principle, al- 
lows to identify the system's different phases and to an- 
alytically predict the phase diagram. We summarize our 



main findings in Section VI and provide a brief conclu- 
sion. 



II. MODEL 

We examine a driven diffusive system which serves as 
a minimal model for the transport on two parallel lanes 
which are coupled by a repulsive short-range interaction. 
The same model describes classical driven spin transport 
with Coulomb blockade. In the following we specify the 
dynamics in detail, presenting both the two-lane and the 
spin-transport representation. 



values or 1, corresponding to a vacant or an occupied 
site, respectively. 

Particles enter from two entangled reservoirs located 
at the left-hand side of the system. At each time step, 
the reservoir is in one of the four possible states: (i) dou- 
ble occupation with probability n* , (ii) only the upper 
reservoir is occupied with probability p| cs — k* , (iii) only 
the lower reservoir is occupied with probability p^ cs — k* 
and (iv) both reservoir sites are empty with probability 
1 — Pres — Prcs + K * ■ Thus, pj es and p^ es are the average 
densities on the upper and lower reservoir, respectively, 
and k* corresponds to the double-occupation density in 
the reservoir. 

In bulk, particles move unidirectionally to the right. 
Due to obstruction, the hopping rate thereby depends 
on the particle configuration at the other lane. A parti- 
cle attempting to proceed by one site is obstructed if a 
particle resides on the subsequent site of the other chan- 
nel. However, this obstruction is relevant only when the 
particle does not experience obstruction in its current po- 
sition, meaning when its current neighboring site on the 
other channel is empty. We model this effect by reducing 
the hopping rate to a value e~ u (U > 0) in this case, 
while the rate is unity for all other configurations. In the 
biological context of molecular motors walking along mi- 
crotubuli, the reduced hopping rate corresponds to the 
spatial obstruction stemming from large cargos attached 
to motor proteins. 

Last, the rules of the model are completed by speci- 
fying how particles leave the system after traversing the 
bulk. Here, we consider that particles at the right bound- 
ary leave with the exiting rates pV, (3^ in contrast to pp. 



A. Two-channel Representation 



B. Spin- Transport Representation 
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FIG. 1: (color online) Illustration of the two-lane representa- 
tion. Particles enter from two reservoirs at the left boundary 
with densities pj ea , pj cs . In bulk, hopping rates depend on 
the particle configuration of the, respective, other lane. At 
the right boundary particles leave at rates /? , /3 . 

Consider particle transport along two parallel chan- 
nels, each of them containing N discrete lattice sites, see 
Fig. [I] Each site may contain at most one particle (on- 
site exclusion), such that the occupation number of site 
i on the upper (lower) channel, n\ (nj), can only take 
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FIG. 2: (color online) Illustration of an exclusion model with 
two internal states, adopting the language of spin transport. 
Particles in state t ii) enter the system at the left boundary 
from a reservoir with density pj cs . (/Ores.) an d leave at the right 
boundary with rate /r (/3" 1 "). In bulk particles hop to the right 
always respecting Pauli's exclusion principle. If an unpaired 
particle moves to a site which is already occupied by a particle 
of the other spin state the hopping rate is decreased to e~ u , 
otherwise, the hopping rate is set to 1. 

The model can be readily interpreted in the context of 
spin transport where it serves as a description for classi- 
cal spin currents. The analogy to the two-channel picture 
is the following: a particle situated at the upper (lower) 
lane maps to a particle with spin up (spin down) . At the 
left boundary particles enter from a spin reservoir with 
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densities pj es , py es . Having traversed the lattice, they 
leave the system at the right boundary with exiting rates 
/Tr, p*. In bulk the particles move to the right always 
respecting Pauli's exclusion principle, i.e. only one parti- 
cle per internal state is permitted per site. According to 
the two-lane representation, the hopping rates depend on 
the particle configuration of the system. A short-range 
repulsive interaction reduces the hopping rate for an un- 
paired particle onto a site which is already occupied by a 
particle of the other spin state to e~ u , as compared to 1 
for the other configurations. The parameter U > may 
be viewed as an effective interaction potential, originat- 
ing from a repulsive Coulomb interaction, where particles 
on the same site (though different spin states) gain po- 
tential energy. In this context, one can also consider an 
increased hopping rate away from a double occupation, 
yet, one can show that this does not change our results 
qualitatively |40j . 

For clarity, we employ only the two-channel represen- 
tation in the following. 

A similar model as the one introduced above was re- 
cently proposed in pQ. However, only symmetric situa- 
tions where considered with entrance/exit reservoirs that 
were equal an both lanes. As a further difference to 
our model the authors did model the exiting processes 
through reservoirs at the right side instead of exiting 
rates. Because we explicitly investigate the asymmet- 
ric case, with entrance/exit properties that differ for the 
two lanes, and because of our usage of exit rates we find a 
multitude of new phenomena, summarized in Section [VI] 
The asymmetry between the two lanes requires a two- 
dimensional generalization of the extremal-current prin- 
ciple. The derivation of this two-dimensional extremal- 
current principle constitutes a key result of our work; we 
show how it successfully describes much of the system's 
behavior. 



III. CLASSIFICATION OF THE SYSTEMS 
SENSITIVITY ON THE POTENTIAL 

The steady-state bulk densities p\ — (nj), p\ = (nj), 
where (•) indicates a coarse-grained time average, con- 
stitute key observables. Because of particle conservation 
their temporal evolution can be obtained from the parti- 
cle flux ji-i onto site i and the one away from it, jf. 
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The currents, jj , jf, contain correlations between neigh- 
boring sites on the lattice. To find an analytic description 
these correlations have to be accounted for by a suitable 
closure relation, e.g. by a mean-field approximation or a 
one-site cluster approximation. 

Stochastic simulations provide another route to gain 
insight into the system's behavior. In this Section we first 
detail the simulation algorithm, and then describe three 



classes of behavior that emerge for different interaction 
strength. 



A. Stochastic Simulations 

We have determined the system's stationary state via 
stochastic simulations with random sequential updating, 
using the dynamic rules introduced in the previous sec- 
tion and employing the Gillespie algorithm [1TJ U2] • We 
have performed time averages over about 10 5 time inter- 
vals, each containing 10 x L time steps and the lattice size 
is set to L = 1000. At the left boundary, the reservoir 
dynamics is specified in terms of the three parameters, 
pj. cs , p^ cs and re*. Here, we restrict the discussions to 
the case of relaxed reservoirs, where the correlations in 
the reservoirs reflect the ones in bulk, which is particu- 
larly illuminating and amenable to a theoretical descrip- 
tion. Then, the double occupation density can be deter- 
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Eq. ( A4 1 derived in the Appendix. In general, and apart 



mined from the average densities, pj. cs , pj cs according to 



from boundary effects such as boundary layers, we found 
constant density profiles in the system. To determine the 
corresponding value of the average density in bulk for 
constant density profiles, we only considered the 0.2 x L 
sites in the center of both channels. Our simulations con- 
firm to a large extent the analytic approximations which 
are to be discussed in the following sections. 



B. Dependence on the Interaction Strength U 

In the case of vanishing coupling, i.e. U = 0, the 
system simply corresponds to two uncoupled TASEPs. 
In the presence of obstruction, and thereby coupling 
between both channels, this picture changes drastically 
upon increasing the effective interaction strength. Our 
stochastic simulations show three regimes of qualitatively 
different behavior, which are illustrated in Fig. [3j In the 
following, we discuss these regimes, and provide a map- 
ping on TASEP for two of them. 



1. Weak Coupling 

In the first regime (I), for small coupling strength U, 
the system almost behaves like two uncoupled lanes, ex- 
cept for the fact that the densities are slightly reduced. 
This regime can be well described by a simple mean- 
field approximation, where correlations between differ- 
ent lattice sites, i ^ j, are neglected, (nj^rat^) rj 

(n? or by a one-site cluster approximation 

which is to be discussed in detail in the following sec- 
tion. In this regime the phase behavior qualitatively cor- 
responds to the one already known from TASEP. On a 
quantitative level, differences arise as the phase transi- 
tion lines are shifted compared to the uncoupled case. 
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FIG. 3: (color online) Average bulk densities obtained by 
stochastic simulations on the upper (red, grey) and lower lane 
(black) for varying potential strength U. The parameters are 
pJcs = Pros = 0.5, ft = 0.1 and ft = 0.3. Three regimes of 
qualitatively different behavior emerge. In the first one (I) 
the system qualitatively behaves like two uncoupled systems. 
In regime II, the density in the lower channel strongly de- 
creases, while the density in the upper channel is still large. 
This nontrivial behavior is discussed in detail in Section IIVI 
and Section [V] In regime III, for strong coupling, the system 
behaves like a one-channel TASEP. 



3. Strong Coupling 

A further increase in U leads to strong obstruction 
between the lanes with new qualitative behavior. In this 
regime of strong coupling (III) double occupancy of a site 
almost never occurs due to the vanishing hopping rate 
e~ u that would yield this configuration. Therefore, the 
transport properties of two coupled channels is similar 
to a single-channel TASEP. This mapping cannot always 
be performed if the boundaries are reservoirs. It then 
only holds for the special case of small reservoir densities. 
Otherwise correlations between the upper and the lower 
reservoir are large. They are transported in the bulk and 
there destroy the effective one-lane behavior. Especially, 
the behavior shown in Fig. [3] would drastically change, 
if the exiting rates ft — 0.1 and ft = 0.3 would be 
replaced by the corresponding right reservoirs, namely 
the description holding for regime II would then apply 
and the bulk densities would not show the one-channel 
behavior. 

Introducing the total density, tj = (nj + n\), and per- 



forming the limit e 



0, which yields 
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the expression for the currents (Eq. Al I, one can identify 
the following mean-field current-density relation, already 
familiar from TASEP, 



J = t{X-t). 



(2) 



2. Intermediate Coupling 



In regime II, the one of intermediate coupling strength, 
an intriguing new phase behavior emerges. For instance, 
for the set of parameters shown in Fig. [3j the density 
in the upper channel remains rather undisturbed by the 
obstruction, while the one in the lower channel drops to 
a comparatively small value. The nontrivial behavior in 
regime II is caused by the influence of the potential on 
the transport properties in bulk as well as on the bound- 
aries. Especially, the exiting current is strongly influ- 
enced by the interaction potential resulting in smaller 
densities at the right boundary than expected for uncou- 
pled systems. We rationalize this behavior in the fol- 
lowing section. Because the systems operates far from 
equilibrium, this change in the boundary conditions has 
a strong impact on the system. 

Secondly, the transport properties in bulk also react 
sensitively to the coupling. As discussed in detail in the 
following section, this causes changes in the currents' de- 
pendences on the bulk densities for potentials larger than 
a critical strength Uc — In 4 ss 1.4. 

Furthermore, other intriguing phenomena, such as do- 
main walls between a low and a high density phase, are 
found. Also, there exist phases where the system depends 
sensitively on both boundaries, namely, the total current 
of the system is fixed by the right boundary while the 
exact value of the densities in bulk depends on the left 
boundary. 



Due to the large potential in this regime, particles on dif- 
ferent lanes are not able to 'overtake' each other. Hence, 
the ratio of the densities in both channels is fixed to the 
value given by the reservoir densities at the left boundary, 
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In Fig. [3] we have used equal reservoir densities for the 
upper and the lower lane, and as a consequence the bulk 
densities of both lanes are equal. 

The exact phase behavior can be determined by relat- 
ing the boundary conditions of the two-lane system to 
the corresponding boundary conditions of the effective 
one-lane TASEP, for which the exact phase diagram is 
known. For reservoir densities p| cs , pj CB < 0.5, the ef- 
fective entering rate is obtained by simply adding both 
reservoir densities: 
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(4) 



The effective exiting rate displays a more complex depen- 
dence on the boundary processes because the individual 
exiting rates influence the exiting current on both chan- 
nels. To find a good estimate of the effective exiting rate, 
we consider the average time a particle spends on the last 
lattice site before it leaves the channel. This time is the 
inverse of the corresponding exiting rate. The weight of 
both waiting times is given by the ratio of particles in the 
upper and lower lane. Hence, a fraction yoJ es / (pj es + p]: es ) 
of all particles spend 1/ft time units on the last lattice 



5 



site, and a fraction pl es / (pj cs + pl cs ) of the particles 1 /(3^ 
time units. The average time is the sum of both times 
weighted with their frequency, 
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yielding the effective exiting rate, 



/3eff = 
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For reservoir densities pj a 



Pics > 0.5, 



the double oc- 



cupation density at the reservoir does not vanish and is 
transported into the system. Hence, the system can ex- 
hibit total densities larger than one if its bulk behavior 
is determined by the left boundary. In this case, the 
description we introduce below for the regime of inter- 
mediate coupling applies. 



where pW)(l — ptw) is the particle current known from 
TASEP and p — k — p^ p^ is the correlation correction 
reducing the current compared to the case without any 
coupling. Here, the double occupation density k is the 
positive solution of the quadratic equation, 



= (1 - e- u )n 2 + [1 - (1 



)p] « 



(7) 



(5) where p = p^ + p^ is the total density. We will employ 
these results in the following section within the frame- 
work of an extremal-current principle to investigate the 
phase behavior as a function of the coupling strength. 

In agreement with the considerations in the previous 
section, the double occupation density vanishes in the 
limit U — > oo , while it simplifies to K = p^ p^ for U — > 0. 
The latter limit connects the one-site cluster approxi- 
mation to the simple mean-field approximation which is 
accurate for small potentials. 



IV. CURRENT-DENSITY RELATION 



The interaction between neighboring particles directly 
affects the transport properties of the system. The cur- 
rent's dependence on the bulk densities is very sensitive 
on the coupling. Above a certain value the current- 
density relation changes qualitatively resulting in new 
phases and phase transitions as we show in the following 
section. 



A. One-Site Cluster Approximation 

With increasing coupling strength U, the occupation 
numbers of the same site on different lanes become more 
and more correlated and a simple mean-field approxima- 
tion fails. However, by employing a one-site cluster ap- 
proximation we obtain a valuable expression for the cur- 
rents depending on the bulk densities as demonstrated 
in PP- To account for correlations between the same site 
on different lanes, we introduce, besides the single parti- 



cle densities p\ 



Pi 



the double occupation 



density on site i, m :— (n[n\) as a new variable 03] . 

Then, the probabilities for the other three particle con- 
figurations on site i, unoccupied or occupied by one parti- 
cle either on the upper or the lower lane, can be expressed 
in terms of p], p\ and k%. We neglect all other correla- 
tions and employ the standard decoupling approximation 
scheme there. 

Assuming spatially homogeneous density profiles, the 
currents on the upper and the lower lane can be expressed 
in terms of p' , p^ and k. The details of the calculations 
are presented in Appendix [XJ We obtain, 



/ = p T (1 - p T ) + M , 
j l = P l (1 - P l ) + P , 



B. Dependence of the Currents on the Potential U 

The current displays a sensitive dependence on the 
strength of interaction, U. For small interaction strength, 
the currents on each lane are almost independent of the 
density on the other lane. Indeed, these currents are ap- 
proximately parabolic with respect to the density in the 
respective lane, see Fig. 4(a) The maximal current on 
the upper channel occurs for p' = 1/2, p^ = and 
p^ — 1/2, p^ = 1 since the obstruction does not affect the 
transport on the upper lane for these densities. With an 
increase in the potential U, the particle flux decreases, in 
particular, for densities around p^ = p^ = 1/2, as shown 
in Fig. M 

This behavior is also reflected in the total current as 
shown in Figs. 4(c)| and |4(d)| For small potentials, the 
total current in bulk displa ys a s ingle maximum located 
1/2, see Fig. |4(c) 



at p> 



In this regime the po- 



(6) 



tential only affects the value of the maximum but not 
its position, i.e. it does not change the topology of the 
phase diagrams. Beyond a critical value of the potential, 
Uc = In 4 « 1.4, the total current displays a qualita- 
tively different behavior, as has been described in Ref. pQ. 
At the critical value, two maxima separated by a sad- 
dle evolve in the current-density relation for the total 
current, see Fig. |4(d) The location of these maxima 
is evaluated in App. |A| Upon further increasing U, the 
maxima move apart. In the limit of large potentials two 
elongated maxima evolve and the saddle becomes a valley 
located at p^ = 1 — p^. The bimodal structure leads to a 
richer phase diagram than in the weak coupling regime, 
U < In 4. Similar extrcma in the currents were found pre- 
viously for one-channel systems, e.g., when next-nearest 
neighbor interactions are taken into account [44| . 
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(a) Single-channel current, y , (b) Single-channel current, y , 
for U = 0.3 for U = 2 




(c) Total current, J, for U = 0.3 (d) Total current, J, for U = 2 

FIG. 4: (color online) Contour plots for the individual (top) 
and the total current (bottom), depending on the bulk den- 
sities, for U = 0.3 (left) and U = 2 (right) using Eq. On 
grey scale black corresponds to vanishing current, and white 
to the respective maximal currents j max , J max . The colored 
(grey) contour lines indicate currents of value 0.95 ■ jmax an d 
0.95- J max (red), 0.8-j max and 0.8- J max (blue), 0.5-j max , and 
0.5- J m ax(green), 0.3-j max and 0.3- J max (orange) from inside to 
outside. Increasing the obstruction strength, the single max- 
imum splits into two, separated by a saddle. The transition 
happens at a critical value of the interaction, Uc = In 4 ~ 1.4. 



C. Influence of the Potential U on the Right 
Boundary 




FIG. 5: The dependence of the bulk densities, p^ = p^ , on 
equal exiting rates /3' = 0* for systems determined by the 
right boundary. It is obtained by evaluating Eq. (JSJ. The 
interaction strength U — 3.0 is above the critical value Uc- 

cles increases. In particular, for potentials larger than Uc 
the double-maxima structure of the bulk current comes 
into play and causes a discontinuous dependence of the 
bulk density on the exiting rates. Such a jump in the 
densities is exemplified in Fig. [5] for the case of equal 
exiting rates. 



V. PHASE BEHAVIOR 

The system's phases depend on the boundary condi- 
tions as well as on the strength of internal obstruction. 
In the following, we employ the extremal-current princi- 
ple to evaluate the analytic expressions for the currents 
obtained in the previous section to get insight into the 
phase behavior. We then discuss these phases and point 
out some special features arising from the coupling. 



At the right boundary (i = L) exiting rates control the 
currents out of the system: 

Jlx = PpI 

Jex = PpI (8) 

These currents are either determining the system or are 
virtual currents which are important for predicting the 
phase behavior in the system as explained in the follow- 
ing section. If these exiting currents also set the bulk 
currents, i.e. if — jg X and = j^ x , we can com- 
pute the bulk densities (which then equal the densities 
at the right boundary) depending on the exiting rates 
via Eq. The densities at the right boundary also 

play a key role for determining the phase diagrams as we 
will see below. For a small coupling strength, we find 
ptU) « 1 — as familiar from TASEP. The transport 
properties change rapidly when repulsion between parti- 



A. Extremal-Current Principle 

The extremal-current principle (ECP) often governs 
the phase behavior of driven diffusive systems [13l HU 
H5] . So far the ECP has only been established for 
one-dimensional systems. Here we describe a two- 
dimensional generalization and show that it successfully 
describes the phase behavior of our two-lane model. 
We start with a short description of the standard one- 
dimensional ECP, and then extend it to two coupled 
lanes. 

The ECP for transport on a single lane can be formu- 
lated by considering two characteristic velocities. The 
first is the collective velocity, v c — dj/dp, which re- 
veals information of the stability of a given bulk den- 
sity p against perturbations: only densities with v c = 
(as well as those determined by the boundaries) are sta- 
ble. The second quantity is the shock velocity v s = 
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\j(pi) ~ HP2] I {pi — p-i) that gives the direction in which 
a domain wall between two densities, p\ and p2, travels. 
In this way, v s determines which of both densities, p\ or 
P2, dominates. To find the system's bulk density, it there- 
fore suffices to first identify the stable densities, using the 
collective velocity, and then, by pairwise comparison via 
the shock velocity, single out the bulk density. These 
considerations are summarized by the extremal-current 
principle: 

j = max j(p) for p+> P-, 

p£[p+,P-\ 

j = min j(p) for p+ < p-, (9) 

where p + is the density at the left boundary, and p_ is 
the density at the right boundary. Hence, the system 
is either determined by the entering or exiting current 
or by an extremal current corresponding to a density in 
between the boundary densities. 

On two coupled lanes, the currents in bulk are gener- 
ically influenced by both lanes. We therefore have to 
consider the dependence of the currents on both p' and 
p^ . As in the one-dimensional case, either the maximal 
or the minimal (total) current (see Fig. [6] blue and green 
area) determines the transport in the system, and the 
velocities v c and v s govern which of both scenarios is re- 
alized. However, in order to decide which of both cases 
applies it is no longer sufficient to compare the densities 
at the boundaries. Because of many potentially conflict- 
ing cases a rigorous derivation of the ECP provides a 
considerable challenge. We have, however, observed that 
the following intuitive version of the ECP describes our 
model's phase behavior in the full parameter space. 

In the one-dimensional ECP an extremal current be- 
longing to a density in the interval determines the sys- 
tem. In the two-dimensional scenario the interval is re- 
placed by an rectangle bounded by the boundary densi- 
ties, {p\-,Pi) and (pJcs, Pro S ). Depending on the bound- 
ary conditions either the minimal or the maximal current 
gives the bulk currents. The currents, which have to be 
considered are the entering and the exiting current or a 
mixture of both, i.e. one lane is determined by the left 
while the other one is determined by the right bound- 
ary. The exiting current can be calculated by equat- 
ing Eqs. ^ and Q with |7]). For the specific example 
[7 = 3, pT, /3^, the corresponding contour line of the 
current-density relation consists of several disjoint lines, 
see Fig. [6] where the red line has four parts which are 
denoted as 1,2,3 and 4. Then, the part that includes the 
point (,0^,/}^) marks the boundary where the maximal 
or minimal current is selected. In the example discussed 
here this boundary is given by line 2 in Fig. [6j The ex- 
trema to be considered for the ECP are located at the 
boundary of the rectangle given by {p[,p l L ) (p$ es) pj cs ), 
or at the extremum on the rectangle. This extremum can 
be inside the rectangle or also at its boundary. 

To illustrate the extremal-current principle we consider 
a path where the reservoir densities are gradually in- 
creased along the diagonal p rcs . = pl cs — Pros for fixed 




FIG. 6: (color online) Contour plot of the total current, de- 
pending on the bulk densities. The bulk densities emerging 
for special values of equal reservoir densities (pl cs = pj cs ), 
namely continuously increasing from to 1, are displayed as 
black dots. The other parameters are 0^ — — 0.3 and 
U = 3.0. The red contour line corresponds to the total exit- 
ing current as emerges if the right boundary determines the 
bulk densities. It marks the transition from the minimal (blue 
area, lower left corner) to the maximal (green area, upper 
right corner) current principle and two phase transitions. In 
the lower left corner the system is in the IN/IN phase, crossing 
the red contour line it enters the EX/EX phase. Upon further 
increase in the reservoir densities the IN/IN phase is reached 
again in the upper right corner, before the MC/MC phase 
arises where the densities are limited by the bulk properties. 

exiting rates and interaction strength. Even though we 
choose this path as an example, our results hold for ar- 
bitrary boundary conditions as exemplified in the follow- 
ing. Fig. [7] displays the bulk densities and the phase 
transitions that occur. For small reservoir densities the 
minimal current is selected, which is given by the left 
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FIG. 7: (color online) Average bulk densities for increasing 
reservoir densities p rea = Pres = Pres- The parameters are set 
to j3' — = 0.3 and U = 3; the situation corresponds to 
Fig. [6] The phases introduced here are discussed in Sec. |VB] 
Red (grey) circles denote simulation results; black lines cor- 
respond to analytical predictions from the ECP. 



reservoirs there. Upon crossing the current contour line 
(1) for the first time, the left reservoir current is larger 
than the exiting current, and transport is determined by 
the exiting rates. Crossing the contour line that includes 
the point {p\,Pl) (^)j does not give rise to a phase tran- 
sition, since now the maximum current determines the 
bulk current. In this domain the exiting current is larger 
than the left reservoir current. The next phase transition 
happens when the current at the left boundary exceeds 
the exiting current. Again, the phase transition occurs 
at a segment of the red contour line (line 3) . Upon a fur- 
ther increase, the second maximum of the current-density 
relation is reached, and the maximal-current phase is en- 



tered. For reservoir densities larger th an ( p\,pY 



i.e. 



for the maximum bulk current, see Eq. (A6) in the Ap- 
pendix, a maximum current is attained. 



B. Phases 

As discussed above the system either adopts its mini- 
mal or maximal current, depending on the boundary con- 
ditions. These extremal currents can be either given by 
one boundary or by an extremum of the current-density 
relation itself. Hence, we can distinguish two classes 
of phases in the system, boundary- and bulk-induced 
phases. The first one is highly sensitive to small changes 
in the boundary conditions, while in the latter one the 
densities arc determined by the bulk properties and do 
not depend on the entering and exiting parameters. 



1. Boundary-induced phases 

The boundary-induced phases depend either on the 
entering or exiting processes, and we consequently dif- 
ferentiate between IN and EX phases. In our model, we 
employ particle reservoirs at the left boundary, but ex- 
iting rates at the right one. As a consequence the left 
and the right boundary influence the bulk densities in 
qualitatively different ways. Indeed, in the case where 
both lanes are determined be the left boundary (IN/IN 
phase), the bulk densities are given by the reservoirs den- 



sities, p 



t 



pL, P l 



pL 



In contrast, if a system is in 



the EX/EX phase, only the total current is fixed to the 
value given by the right boundary, whereas the bulk den- 
sities also depend on the reservoir densities at the left 
boundary, for p[ cs pj cs holds ^ p^ . Further mixed 
phases (IN/EX or EX/IN) may become relevant where 
the current on one lane is determined by the left bound- 
ary and the current on the other is fixed by the exiting 
current. 

The IN/IN and the EX/IN phase are exemplified in 
Fig. [8] where we show the average bulk densities depend- 
ing on the upper reservoir density p| os . Parameters are 
chosen such that a transition from the IN/IN to the 
EX/IN phase emerges at a certain value of pj cs . Only in 




FIG. 8: (color online) Average bulk densities for different up- 
per reservoir densities pi es .. Parameters are pj cs = 0.02, /3' = 
I3 l = 0.3, U = 3.0. If both reservoir densities are small, the 
system is in the IN/IN phase where both bulk densities are 
given by the reservoirs. This changes for larger pl cs ,. While 
the lower lane is still determined by the left boundary, the 
upper lane is governed by the right boundary (EX/IN phase). 



the IN/IN phase do the bulk densities vary when chang- 
ing the upper reservoir density pj es ; in the EX/IN phase 
they are almost undisturbed by changes in pj es . 



2. Bulk-induced phases 

In the bulk-determined phases, the current is given by 
an extremum of the current-density relation. Here we re- 
strict the discussion to the case where a maximum occurs 
(MC phase), yet one can find parameter regions where a 
saddle fixes the phase behavior [40 . 

In the maximum-current phase a localized domain wall 
can emerge, separating a high-density regime at the left 
and a low-density regime at the right, see Fig.[9j This sce- 
nario, previously found in Ref. I13|. occurs if both max- 

and (4 / >/>//) 



ima {p\,p\) and (pj T ,p\ T ), Eq. (|A6|) in the Appendix, 



are accessible. Each domain along the lane then corre- 
sponds to one of the maximal currents. In particular the 
current is continuous at the domain wall. In Fig. [9j the 
reservoir density at the left boundary is larger than the 
density corresponding to the second maximum, while the 
density at the right boundary is smaller than the density 
corresponding to the first maximum. Hence, in the vicin- 
ity of each boundary the density which is closer to the 
respective boundary density arises. Depending on the 
exact values of the boundary densities, the domain wall 
is located between the left hand side and the middle of 
the system. Increasing the reservoir densities the domain 
wall moves further to the right. In Fig. [7) the MC/MC 
phase is shown for symmetric boundary conditions. The 
black line denotes the analytically obtained density cor- 
responding to the second maximum of the current density 
relation. It is in excellent agreement with the simulation 
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FIG. 9: Density profile of a state exhibiting maximal current. 
The parameters are pj cs = pj cs = 1, 0^ = fi^ = 1 % U = 
3.0, and L — 1000 . The current remains spatially constant 
while the density shows a high value at the left and a low 
value at the right. Both densities correspond to maxima in 
the current-density relation, pi, pn. An unusual localized 
domain wall forms between them, see text. 



results. 



C. Phase Diagrams 

Employing the extremal-current principle, the phase 
diagrams can be characterized completely also for un- 
equal boundary densities.. As described above, there 
exist several phases which are either determined by the 
boundaries or by the bulk transport properties of the 
system. According to the ECP the phase transition lines 
are given by equating the entering and exiting currents 
or by the structure of the current-density relation. The 
bulk current dominates if the maximum on the rectangle 
bounded by the corners (pl cs , p]: cs ), (p\, p L ) is not given 
by either of these points. It is clear that a transition from 
a left-reservoir or right-exiting-currents dominated phase 
to a maximum-current phase can emerge only if the ex- 
ternum traverses the boundary of the rectangle consid- 
ered. Further transitions to an IN/EX and EX/IN phase 
are identified using again the extremal-current principle 
employing contour lines corresponding to the total cur- 
rent in the mixed phases. 

In Fig. [lO] we exemplify a phase diagram depending on 
unequal reservoir densities. The exiting rates are fixed 
to (ft , /r - = 0.3 and the potential is set to U = 3, a value 
where the current-density relation shows two maxima. 
The parameters are identical to the ones of Figs. [6] and 
[7] In Fig. 10 the stochastic simulations (red dots) are 
in good agreement with the analytic calculations (black 
lines). Because the latter ones were obtained employing 
the ECP also for unequal reservoir densities, the strength 
and generality of the ECP can be confirmed. Due to 
the symmetry between both lanes, the phase diagram is 
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FIG. 10: (color online) A generic phase diagram with the 
reservoir densities as control parameters. The exiting rates 
are fixed to ft' = ^ = 0.3 and the interaction strength is set 
to U = 3.0. The red (grey) dots denote simulation results of 
the transition lines, while the black lines are calculated em- 
ploying the ECP. Besides the combinations of phases which 
are determined by entering or exiting currents, phases which 
do not exist for U < In 4 are also present. Namely a second 
IN/IN phase, which exhibits a high density, and the MC/MC, 
where the maximal current determines the bulk densities, oc- 
cur. 



symmetric along its diagonal, pl es = p^ es . 

In the lower left corner, where both reservoir densities 
are small, both channels are determined by the enter- 
ing current (IN/IN phase). Increasing only one reservoir 
density, the minimal current is no longer given by the en- 
tering current on the respective lane. Hence the system 
reaches the EX/IN phase (IN/EX phase). The phase- 
transition lines between the IN/IN and IN/EX phase, 
as well as the ones between the IN/EX and the EX/EX 
phase, are given by contour lines of the current on an indi- 
vidual lane, Eq. ^ . In the middle of the phase diagram, 
both channels are determined by the exiting currents. As 
mentioned above, only the total current is fixed in this 
phase, in contrast to the bulk densities which also de- 
pend on the reservoir densities at the left boundary. The 
dotted line marks the region where the bulk densities on 
both channels are equal. On this line the bulk densities 
can be calculated employing Eq. Q. For a further in- 
crease in the reservoir densities, the IN/IN phase arises 
again where the bulk currents as well as the bulk den- 
sities are given by the left boundary. This phase would 
not arise if we had chosen entering rates, rather than par- 
ticle reservoirs, at the left boundary, because boundary 
densities larger than 1/2 would then not emerge. In our 
case, such large densities cause currents larger than the 
entering current. These dominate the system according 
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to the extremal-current principle. For even larger reser- 
voir densities the crest of the current-density relation, 
which marks the phase transition from the IN/IN to the 
MC/MC phase, is reached, and the bulk densities remain 
constant at the values where the second maximum is lo- 
cated. 



VI. SUMMARY AND CONCLUSION 

In this paper we have examined a driven two-channel 
model where the motion of particles along both chan- 
nels is coupled via a repulsive short-range interac- 
tion. The latter causes intriguing new phenomena and 
phases. Varying the strength of particle obstruction, 
three regimes of qualitatively different behavior evolve. 
First, for small coupling the system approximately acts 
as two uncoupled lanes, i.e. the phases and phase dia- 
gram qualitatively correspond to the ones already known 
from TASEP. These results can be confirmed by employ- 
ing a simple mean-field approach or a one-site cluster ap- 
proximation. Second, for moderate interaction strength 
the transport properties of the system are strongly in- 
fluenced by the obstruction between neighboring parti- 
cles. This regime emerges around a critical interaction 
strength, Uc = In 4 pQ, where a single maximum of the 
current-density relation splits into two, separated by a 
saddle. Third, when the obstruction is large, the two 
coupled lanes effectively behave as a single one. In this 
case, we have identified a mapping from the parameters 
governing entering and exiting processes in our system 
to effective rates for a single-lane TASEP. This mapping 
then allows to carry over known results from TASEP such 
as its phase diagram. Hence, for different strengths of the 
obstruction a variety of peculiar phases surface which can 
be accessed by manipulating the system at the bound- 
ary only. The boundary-sensitive phases respond grad- 
ually upon tuning the left reservoir or the exiting rates, 
whereas the maximal-current phase is robust against such 
changes. 

In contrast to previous work pQ, we explicitly investi- 
gated the two-channel system with a broken lane symme- 
try. We thereby followed the proposition in pQ that the 
ECP might be generic for multichannel system. As a key 
result we indeed derive a suitable generalization of the 
ECP to two dimensions. We thereby approve that not 
the densities but the currents govern the transitions be- 
tween the minimal and the maximal current principle, a 
distinction that cannot be made within one-dimensional 
or symmetric situations. The accuracy with which the 
system's phase behavior can be predicted with our two- 
dimensional ECP is astonishing and suggests it for fur- 
ther applications. In further contrast to pQ we specifi- 
cally investigated the dependence on the potential U. For 
small potentials the system's behavior is not sensitive on 
whether reservoir or rates are chosen at the boundaries. 
However, for intermediate and large potentials, reservoirs 
or rates at the boundaries make a difference. For exam- 



ple, the system does not behave like a one-channel system 
for large potentials if reservoirs are chosen instead of rates 
at the right boundary, as has been done in pQ. Because 
we employ exiting rates at the right boundary, we observe 
novel effects such as a strong impact of the potential U on 
the exiting current (Fig. [5]), and two transitions instead 
of one in the bulk densities upon increasing U (Fig [3]). 

It would be interesting to consider also lane changes 
(resp. spin flips) of the particles as they proceed along 
the channel. Such events are expected in realistic appli- 
cations such as intracellular transport, highway traffic on 
multiple lanes, or hopping transport of electrons through 
a chain of quantum dots. The correlations induced by fre- 
quent lane changes [3H |3S] , ignoring mutual obstruction, 
are accurately described by a one-site cluster approxima- 
tion similar to the one we employed. Yet, for rare lane 
changes a simple mean-field approximation suffices to de- 
scribe the arising localized domain walls [T^]. The com- 
bination of obstruction, lane switching and possibly also 
defects constitutes a promising route to discover novel 
and unexpected collective phenomena in driven trans- 
port. 

In conclusion, we have shown that the ECP can be 
generalized to higher dimensions to serve as an appro- 
priate tool for the investigation of driven multi-channel 
systems. We therefore think that extremal-current prin- 
ciple is a promising starting point to achieve a deeper 
understanding of complex transport phenomena. 
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Appendix A: Current-Density Relation 

In this Appendix, we derive an analytic description 
for the current-density relation within a one-site cluster 
approximation. 

The individual currents onto a site i can be obtained 
be evaluating the particle fluxes onto this site, 

31= (e~ U ("Li ~ nUnU) (4 - 44) ) 

+ ( - 4-i4-i) (i - 4 - 4 + 44) ) 
+ (4-i4-i(i-4)), (ai) 

and similarly for jf. This expression is evaluated em- 
ploying the one-site cluster approximation discussed in 
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Sec. IV A The essence of the approach consists of con- 



sidering all four possible states of the two opposing sites, 
whereas all correlations between neighboring sites are fac- 
torizcd. Thus, a complete description is achieved in terms 
the mean double-occupation density ft, = (n[n\) besides 
the average densities p\ = (nt) and p^ = (nj). Then the 
closure relation for the currents is derived 



+ {p\-i ~ ^-^{l - p\ - p\ + Hi) 

+ Ki_i(l-pt). 



-t 



(A2) 



Similar to Eq. ( Al ), the time evolution for the double- 



occupation density Ki involves averages of products of 
occupation variables n,. Within the same truncation of 
the hierarchy, one finds, 



d t Ki 



-u ( J 



+ e~ (pi_ t - «i_i)(p! - Ki) 

+ Ki-l(Pi+Pi ~ 2«i) 

- «i(2 - pLi - pLi) • 



(A3) 



In the steady state and for spatially homogeneous den- 
sity profiles, the double occupation density k obeys a 
quadratic equation with solution, 



Xp-1+ y/(\p- I) 2 +4A(1 - A)pV 



2A 



(A4) 



Here, A = 1 



corresponds to the inverse dwell time in 



the obstructed case and p ■■ 



p^ to the total density. 



Combining this result with the closure relation for the 
current, Eq. (A2), the closed expression for the current- 



density relation follows, 



f = p t {l-p' t - p X ) + K, 
j X = p^ (1 - p^ — p X ) + K . 



(A5) 



The total current J = + displays a single a max- 



imum for small interactions U located at p^ = p^ 



1/2- 



This maximum is replaced a saddle for strong coupling 
and two maxima of equal height appear on the diagonal. 
These new maxima are located at, 



Pi,u = Piji = 2 ± 4(1 _ e -U) • ( A6 ) 



These solutions are only real for potentials larger than 
the critical value, Uc = hi 4. 
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